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INTRODUCTION 


This report covers Mork done on NASA Contract NAS8-33691, "Investi- 
gation of Electrodynamlc Stabilization and Control of Long Orbiting Tethers." 
In addition to the original topic of studying how electric currents could 
be used to control the motion of a long wire In space, the report also In- 
cludes some new techniques for controlling Instabilities during retrieval 
and an Investigation of how the tether can be used for launching satellites 
from the Space Shuttle. The appendices to this report present derivations 
and analyses of a general nature used In all of the areas studied. 
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stabilization 


1.1 Introduction to Electrodynamic Stabilization 

The objective of this study is to devise an algorithm for using electric 
currents to control pendular oscillations of the Skyhook wire induced by 
various perturbing forces. The electrodynamic force is perpendicular to the 
wire and the earth's magnetic field. The principle limitation of this techni- 
que is that the sign of the wire current (positive or negative) and the 
direction of the earth's magnetic field cannot be controlled so that the 
method offers less flexibility tnan thrusters. Methods have been devised for 
taking best advantage of favorable geometries for electrodynamic damping of 
tether oscillations. 

1.2 Electrodynamic Equations 

A wire of length dl moving at velocity v with respect to a magnetic field 
B experiences an electromotive force V given by v x B • |dl(. The current I 
generated by this electromotive force causes a force F « I x B |d1| on the 
wire. This force is always a drag force unless the orbit is so high that the 
rotational velocity of the magnetic field is greater than the orbital velocity 
of the wire. The maxiimim current that can flow is V/Z where Z is the total 
impedance of the circuit* including the wire resistance, the contact resistance 
of the electrodes with the plasma, and the effective impedance of the external 
circuit in the plasma. The damping algorithm assumes that the desired current 
is within the maximum possible current. In situations where it is not, the 
effectiveness of the dancing would be reduced. 

1.3 Transverse and Vertical Forces on the Tether 

The problem of controlling tether oscillations arises from the fact that 
the Shuttle has direct control over the tether system only through reeling or 
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unreeling the wire. Apart from accc^npanylng Coriolis forces, this method 
allows control of only the vertical position of the subsatelllte. Coriolis 
forces cause the subsatelllte to move In the direction of the orbital motion 
when the subsatelllte moves closer to the earth, and to lag behind the Shuttle 
when the motion 1$ to higher altitudes. There Is also weak coupling between 
the vertical and out-of-plane movements. However, damping occurs only on 
deployment, and reeling In the tether causes negative damping. Other forces 
acting transverse to the wire are the restoring force of the gravity gradient, 
and atmospheric drag. One can, of course, have active control of subsatelllte 
motion using small rockets. If the wire Is an electrical conductor, and cur- 
rents are allowed to flow, there will be an electrodynamic force on the wire 
which Is perpendicular to the wire and the direction of the earth's magnetic 
field. The electrodynamic force has the attractive property of allowing trans- 
verse forces to be applied to the wire without the expenditure of fuel or 
energy. 

1 .4 Tether Instability Modes and Causes 

Tether Instabilities can result during rapid retrieval as a result of 
Coriolis forces and residual angular momentum of the tether system with respect 
to the Shuttle. During station keeping, variations In atmospheric drag can 
excite pendular oscillations particularly when there Is resonance between the 
natural frequency of tether oscillation and the frequency of the drag varia- 
tions. The drag variations of principle concern ar« the diurnal bulge caused 
by solar heating of the atmosphere, the equator to pole variation In density 
arising from earth oblateness, and the corotation of the atmosphere as the 
earth rotates. 

In an equatorial orbit, the diurnal bulge can excite In-plane oscillations. 
Since the perturbing force Is at orbital frequency and the natural period for 
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in-plane oscillations 1: ^ times the orbital frequency, the tether oscillation 
cannot stay In phase with the driving force and beats result. In-plane oscilla- 
tions can also result from orbital eccentricity since the drag on the sub- 
satellite Is larger at perigee than at apogee. 

In a polar orbit. In-plane oscillations can result from the diurnal bulge 
which causes a perturbation at orbital frequency, and the equator to pole varia- 
tion from earth oblateness which causes a perturbation at twice the orbital 
frequency. Beat phenomena would result since neither perturbation Is at the 
natural frequency for In-plane oscillations. 

Gut-of-plane drag forces will act on a subsatellite In a polar orbit as a 
result of the corotation of the earth's atmosphere. The force occurs twice per 
orbit In an easterly direction at the equator crossings. Since the out-of-plane 
natural oscillation frequency Is twice the orbital frequency, this case must be 
examined In detail because of the possibility of resonant excitation of oscilla- 
tions. To do this, let us use an x-y-z coordinate system with the z axis at the 
Northpole, the x axis In the equatorial plane, and the y axis directed to form 
a right handed coordinate system as shown In Figure 1. 

Let us start the Shuttle and subsatellite on the x-axIs with a velocity 
In the z direction. Since the transverse force due to corotation of the atmo- 
sphere Is 1n the direction, let us assume the subsatellite has a pendulum 
velocity In the +y direction In addition to the orbital velocity In the z 
direction. Let the angular position e of the system along the orbit be 
measured north frcxn the x axis. Table 1 lists the out-of-plane position and 
velocity of the subsatellite In the y direction and the magnitude of the 
transverse drag force for the first half orbit. At e • 0, the out-of-plane 
velocity y Is In phase with the perturbing force Fy by design In the Initial 
conditions. At e « 180^ y Is out of phase with Fy. If the magnitude of the 
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Figure 1. Coordinate system for study 
of out*of>p1ane oscillation. 
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Tabic 1. Out-of-planc position and 
velocity of sid)sate111te 
and perturbing force vs«> 
angular position along 
the orbit. 


e 

y 

• 

y 

Fy 

0 

0 



45 

% 

0 


90 

0 

-^M 

0 

135 


0 


180 

0 

% 

-’'h 


corotation force Is the same at the equator crossings on opposite sides of 
the earth, the perturbations oppose and tend to cancel each other. In general 
the forces will not be equal because of the diurnal bulge or orbital eccen- 
tricity. In this case there will be a resonant build-up of the out-of-plane 
oscillation. 

1.5 Daim)1nq of Simple Pendulian Motion 

For small oscillations the horizontal motion of a mass m at the end of 
a rod of length i with a vertical force on the mass and a damping coefficient 

b Is 

•• • 

mx ♦ bx ♦ -J- ■ 0 

The frequency u of the pendulum Is 

-b ± Jb^ - 4fflF /I 

w 


2m 





To obtain critical dmiplng, the quantity under the square root Is set to 
zero to obtain the critical damping coefficient 


b 




The actual Skyhook system Is a spherical pendulum oscillating in a 
rotating frame of reference. If the electrodynamic drag force were In the 
same plane as the pendulian oscillation one could use a damping algorltNn 
similar to the above case for providing an electrodynamic damping force 
during half of the oscillation cycle. The electrodynamic force cannot be 
used over a full cycle unless It Mere possible to reverse the direction of 
the current. In general, the electrodynamic force will not be parallel to 
the velocity and the pendulimi oscillation will not be confined to a plane. 


1.6 Vector Damping 

In Figure 2 an object with velocity v Is acted on by a force F. In 

Figure 2a, the angle e between the velocity vector and the force vector Is less 

than 90**. The velocity v can be resolved Into consents v^ and Vh which are 

perpendicular and parallel respectively to the force F. The force F cannot 

•¥ ■¥ 

change the magnitude of Vj^. Since e Is less than 90**, v^ Is In the same 

direction as F. The force F will Increase the magnitude of v„ and therefore 

Increase the velocity and kinetic energy of the object. In Figure 2b, e Is 

greater than 90* and v^ Is anti 'parallel to F. The force F will decrease the 
♦ 

magnitude of Vy and reduce the kinetic energy of the object. The magnitude of 
Vn can be reduced to zero at which point the angle between v and F Is 90**. 

If F continues to act on v, v will acquire a component of velocity v^ which Is 
parallel to F and the angle between v and F will be less than 90*. The situa* 
tion then reverts to the case In Figure 2a. If the direction and %\gn of F 
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Figurt 2. 


Force F acting rm an object moving with velocity v. 
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art fixtd but tht magnitudt of F can be controlled, the force can be used to 

♦ 

remove energy from the object by setting the magnitude of F proportional to 

the magnitude of as long as v„ Is anti-parallel to F, and setting F to 
■¥ ♦ 

zero when v„ Is in the same direction as F. An algorithm such as this will 

not add energy to the object and will remove energy whenever the geometry Is 

favorable. The discussion above also applies In three dimensions since v can 

♦ 

be resolved Into components parallel and perpendicular to F. Figure 2 can be 

used for three dimensional vectors by choosing axes that lie In the plane 
^ -*• 

determined by v and F. 


1.7 Damping of Spherical Pendulum Motion 

Stabilization of the Skyhook tether requires removing any velocity of the 

tether with respect to the equilibrium velocity of the tether. In Figure 3. 

V 2 Is the velocity of the Shuttle. Vi Is the velocity of the subsatelllte. and 
♦ 

v| Is the velocity of the equilibrium position of the tetiMr. Ttm vector v{ 1$ 
* 

obtained from v^ by the equation 


where r is the distance of the Shuttle from the center of the earth and t 
Is the length of the tether. The velocity to be damped 1$ av given by 

Av ■ Vj - Vj 

The direction of the electrodynamic force on the wire Is 

^ ^ 

F ■ I X Bi 


The vector damping al^rlthm described In the previous section can be applied 
to the vectors F and Av using a suitable damping coefficient b^ as described In 
Section 1.5. The mgnitude of F can th*n be controlled according to the 
equation! 

F ■ -b^ AVu 

♦ ♦ ♦ 
where AVg Is the component of av parallel to F. 



Figure 3. Velocity of the subsatelllte. Shuttle, and 
equillbriun position of subsatelllte. 


0 . 



In general* the electrodynamlc force F and the relative velocity av Mill 
be continuously changing along the orbit of the Shuttle. The geometry will 
be favorable for damping at some times and not at others. The vector damping 
algorltNn Is designed to take the best advantage of opportunities for damping 
on a dynamic basis. 

1.8 Current Generation and Control 

Generation and control of currents In the Skyhook tether can be accomplish- 
ed In different ways. In a coim>1ete1y passive mode electrons In the wire are 
subjected to a force arising frcnn the tether's velocity through the magnetic 
field such that they are driven to the botoom of the wire (If the orbit Is 
easterly). Assuming the wire Is coated, and there are terminating electrodes, 
the electrostatic potential developed at the ends of the wire attracts charged 
particles from the plasma setting up a current loop. The current can be con- 
trolled by Inserting a variable resistance between the wire and the electrode. 

An alternative to a collecting electrode Is a plasma contactor. This would be 
a device that emits electrons or Ions for completing the circuit to the plasma. 
The current could be limited by controlling the emission of electrons or Ions. 
Such a device might be especially useful at the lower end of the wire, which 
collects positive Ions, since collection efficiencies are low for a passive 
negative electrode. If the orbital geometry Is such that the voltage re- 
sulting from orbital motion Is small. It may be possible to drive currents 
with a plasma contactor that accelerates charged particles with a high voltage. 
This could allow electrodynamlc control In cases where the electrodynamlc force 
f * T X ^ I iTj Is non-zero, and In a useful direction, but the driving potential 
f ? Is low either because v x f Is small or the component along t Is 
small. As an exan^le, consider the case shown In Figure 1. If we neglect the 
offset of the magnetic pole, the magnetic field Is parallel to v and therefore 
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there Is no natural driving potential. However If a current could be forced 
In the wire, the electrodynamic force F ■ T x B |I| would be non-zero and 
parallel to the y axis. 

The maximum current that can be obtained depends on many factors, 
particularly In a passive mode. The collection efficiency of electrodes depends 
.on the size and potential of the electrodes and the Ion density In the plasma. 
Collection of electrons Is much more efficient than positive Ion collection 
because of the high thermal velocity of electrons. The resistive limit to the 
current Is the potential difference across the wire divided by Its resistance. 
The potential difference from orbital motion depends on the direction of the 
magnetic field with respect to the orbital veloctly. 

There are various engineering considerations such as the Insulation 
required to withstand the voltages developed In the wire, and methods of pro- 
tecting payloads sensitive to high voltage. Coating the wire Is necessary for 
accurate control of electrodynamic forces since an uncoated wire would allow 
uncontrolled charge collection along the wires. 

1.9 Computer Program for Heuristic Study 

A computer program has been written for numerical Integration of the In- 
plane (e) and out-of-plane U) displacements of the tether vs. time. The 
vector damping algorithm Is Included In the program. 

The program Is designed to run in three modes. In the first mode, there 
Is no electrodynamic force and the tether oscillation results only from atmos- 
pheric drag. In the second mode, the current In the wire Is computed as the 
Induced voltage V divided by the tether resistance R. In the third mode, the 
current I Is determined by the damping algorithm. This last mode assumes that 
the required current can be obtained and Is designed to study what can be done 
In principle with thr electrodynamic force. 
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When the product I x B becomes very small as a result of the magnetic 
field B being nearly parallel to I, the damping algorithm calls for unreasonably 
large currents to compensate. To avoid this, a maximum current parameter was 
added to the program. The damping Is still effective with the current limited 
to a reasonable value since the problem occurs only at certain points In the 
orbit. A feature Is Included In the program to allow the dancing algorithm 
to operate on the vector 




where fg and f^ are weighting factors for the velocity components Vg and v^ 
In the e and ip directions. 

If the damping algorithm Is applied to the out-of-plane component v^ of 
the pendulum velocity of the tether, the total velocity Vj may be either In- 
creased or decreased depending on whether the cosine of the angle ouj. between 
the electrodynamic force F and the velocity Vj, given by 



Is positive or negative. The damping force Is proportional to the cosine of 

the angle between F and v. given by 

♦ 



Some runs were done setting |F| >0. when cos > 0. This reduces the 
opportunities for the damping force to be applied but guarantees that Vj will 
not be Increased. In certain situations this condition on cos causes an 
oscillation to be set up which forces the numerical Integrator to take extremely 
small steps. In effect halting the Integration. The application of the damping 
force causes cos 0 ^ to become more positive. If cos Oj Is passing from a posi- 
tive value through zero because of changes In the direction of the magnetic 
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field B. and a damping force 1$ suddenly applied when cos oj < 0., a feedback 

loop Is created which causes the dwnping force to be turned on and off rapidly. 

The phenomenon Is analogous to the cycling of a thermostat. The problem can 
be corrected either by having a proportional controller, or by Introducing a 
time lag In the feedback loop to prevent rapid oscillation. The basic damping 
algorithm Is a proportional controller since the damping force Is proportional 
to the cosine of the angle between F and the velocity being damped. A time 

lag has been Introduced by setting |F| * 0. for cos ^ and setting |F| 'v 

cos for cos oj <_ Cq||j. By choosing values such as Cqj^ * -.1 and Cgpp » +.1, 
rapid cycling of the damping Is prevented. 

Another problem causing very small Integration steps occurs In near polar 
orbits. The Induced voltage In the tether Is very low and can pass through 
zero slowly as the orbital geometry changes. The damping algorithm was written 
to apply a current which Ignores the value of the voltage but turns on or off 
as the sign of the voltage changes. The Integrator takes very small steps as 
It approaches a change In sign of the voltage. Setting separate voltage cutoff 
values Vgi^ and Vgpp for turning the damping on and off solves the problem. The 
difficulty appears to be one of Integrator accuracy In a grazing Incidence 
approach to a discontinuity. No feedback oscillation was observed. 

1.10 Atmospheric Model 

The atmospheric density p Is given by the equation 

A A 

Pp P • Ps 

0 - pq + -2" 2 x) + pg - — . 

Ip ‘ Pjl 

where pg Is the constant part, Pg Is the equatorial bulge factor, \ Is the 
latitude, pg Is the diurnal bulge factor, p Is a unit vector pointing from the 
center of the earth to the Shuttle, and p Is a unit vector pointing to the 
diurnal bulge. The positions of the diurnal bulge and magnetic pole, and the 
constants pg, p^, and pg are Input parameters to the program. 
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1.11 Magnetic Field Model 


In the case of the earth» we can write the magnetic field B to a fairly 
good approximation In terms of a dipole field centered In the earth's center. 
In reality, the dipole Is displaced from the center (by 300 km). Therefore, 
we can assume that 

B ■ x(yj + zk) + (3x2 . ^ 2)5 

+ — ^ y(zk x 1 ) + (3y2 - r^)J 

r5 rS 

3M_ . * M- 

*jfz(x1 +yj) +-^(3x2 - r 2 )k , 

where 

r2 ■ x2 + y2 -f z2 

X > a cos a 
y « a sin a cos 1 
z » a sin a sin 1 
a > semi -major axis of the orbit 
Also, with reference to Figure 4, 

■ M sin 11® cos X 
My • M sin n® sin x 
- M cos 11® i 0.98M 





Figure 4. Reference system for the 
earth's r/iagnetlc field. 
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1.12 Cases Run In Heuristic Study 


Runs have been made for equatorial, Inclined, and polar or near polar 
orbits. In an equatorial orbit, neglecting the offset of the magnetic pole, 
the tether oscillation and the electrodynamic force F are both In-plane. With 
no damping, the In-plane oscillation due to drag perturbations shoii beat 
phenomena and there Is no build-up of the amplitude with time since there Is 
no resonance. When the damping algorithm Is activated, a current Is applied 
on the forward swing of the tether. This reduces the amplitude of the 
oscillations and keeps the tether closer to the equilibrium position determined 
by the average value of the drag force. Oscillations cannot be completely 
eliminated since there is always a perturbing force present. 

In an Inclined orbit, the drag force and electrodynamic force both have 
In-plane and out-of-plane components. Without dan^)1ng, the In-plane oscillation 
shows beats as expected. The out-of-plane component Is Initially much smaller 
than the In-plane component but grows with time because the drag perturbations 
are resonant with the half orbital period of the out-of-plane oscillations. 

If the damping algorithm Is applied Immediately to the total velocity v of the 
subsatelllte, the In-plane oscillation Is reduced but the out-of- plane oscilla- 
tion Is made larger than It would be initially without damping. If the out-of- 
plane component Is given an amplitude comparable to the In-plane oscillation, 
both the In-plane and out-of-plane oscillations are reduced by applying damping. 
If there Is no Initial oscillation amplitude and the damping algorithm Is 
applied only to the out-of-plane component v^, the out-of-plane oscillation Is 
prevented from building up with time and the dancing currents are small, since 
the driving force In the out-of-plane direction Is small. If the out-of-plane 
amplitude Is Initially large, and dat ing Is applied only to v^, the out-of- 
plane oscillation Is reduced to a small value, but the In-plane amplitude may 
be made larger than It would be without damping. This can be prevented by 
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applying damping to only whan tha angla batwaan F and tha total valocity v 
Is graatar than 90*. Sinca this raducas the opportunities for damping by 
about a factor of two, tha danq>1ng takes about twice as long and the final 
amplitude obtained for the out-of-plane oscillation Is about a factor of two 
larger than that obtained with full time damping. 

In a polar orbit, the electrodynamic force Is In the out-of-plane 
direction If the offset of the magnetic pole Is neglected. The Induced voltage 
In the wire Is small and exists only as result of the rotation of the magnetic 
field with the earth. The vectors I and B are perpendicular at the equator and 
parallel at the poles. If a current can be made to flow In the wire either 
using the naturally Induced voltage In a very low resistance wire, or by using 
voltage sources, the electrodynamic force can be used to control the out-of- 
plane oscillation. This Is significant since It Is the out-of-plane oscilla- 
tion that Is subject to resonance phenomena. 

Monthly Reports 3 and 4 dated December 1979 and January 1980 give the 
details of the runs described above. Figure 5 shows a san^le plot from 
Monthly Report No. 4. The Shuttle Is In a 220 km altitude orbit Inclined 45® 
with respect to the equator. The tether Is 10 km long and deployed upward. 

The out-of-plane angular velocity Is Initially 1.5 x 10"** radlans/sec. and the 
damping algorithm Is applied to the out-of-plane component. Atmospheric drag 
Is not Included. Part (a) of the figure shows the In-plane angle, part (b) 
the out-of-plane angle, and part (c) the current. The first column In each 
graph Is the time In seconds. The orbital period Is 5333 seconds. The second 
column Is the angle (In radians) or current (anps). The out-of-plane oscilla- 
tion Is quickly damped out by the dancing algorithm. 

In all cases described above, the damping algorithm assumes that the 
desired currents can be obtained by seme means. The next section discusses 
practical considerations In generating tether currents. 
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FIGURE CAPTIONS 


Figure 5(a). In-plane angle (radians) vs. time (sec.). 
Figure 5(b). Out-of-plane angle (radians) vs time (sec.). 
Figure 5(c). Damping current (amps) vs. time (sec.). 
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1.13 Limits to tha Ttthtr Current 


In assessing plausible values for the tether current we must take Into 
account, a) the collection of the end electrodes of tN tether (conducting 
surface of the Shuttle at the one end and conducting balloon or subsatelllte 
at the other end); b) the Internal resistance of the tether; c) the Impedance 
of the "circuit" extemil to the tether. 

a) Collection of the end electrodes of the tether 


An analysis of current and potentials In a coated tether configuration Is 
contained In Reference 1 (p. 30 ), TN numbers quoted there have taken Into account 
only the collection properties of the electrodes. All other resistances were 
supposed negligible. Furthermore, a maximum electromotive force aV ■ VBL was 
supposed to be available (aV 20 kv for L ■ 100 km). 

Table 1, Reference 1, contains a sunmary of limiting cases. Here 
by configurations A and B, we mean the cases with the Shuttle upward or 
downward respectively. S stands for Shuttle, b for balloon, the currents 
IjO (J ■ 1. e for Ions and electrons respectively) are thermal currents 
defined by 


1 , 0 . ( 2 ) 
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whert A Is the collecting srea. These are essentially the currents collected 
by an electrode (of appropriate polarity) when Its potential with respect to 
the plasma Is small In comparison with particle thermal energies* I.e .. e^, ^ ^ 
[In the opposite case (ea/KT » 1), the collection efficiency can Increase 
very much]. 

Finally, for the numbers reported in Table 1, the conducting part of the 
Shuttle has been chosen as equivalent to a radius r^ ■ 2.78 m. Looking for 
example at configuration A, L ■ 10^ and with all the available potential on 
the balloon (or subsatelllte) we see that the current 

1 ■ 0.19 amps (3) 

corresponds to the thermal current of ^e Shuttle. 

Increasing the balloon radius does not serve to Increase the current 
unless very large radii are reached; where all the available potential has 
been displaced to the Shuttle* a situation which one clearly wants to avoid. 

Itotlce that the same current value Is also obtained (always In configura- 
tion A) for the tether of L ■ lO^km. The only way to Increase current* with 
respect to the value (3) would be that of using an Ion gun at the Shuttle end, 
either Increase accordingly the balloon dimensions or (In tim case of a sub- 
satellite}* put an electron gun on the subsatellite. 

Let us evaluate the current for a situation where* because of the orbit 
Inclination* the driving e.m. force Is very much decreased. In this limit 
both electrodes can be at very low potentials* so that they are only able to 


draw their thermal currents. In this case we see that, for reasonable 
balloon or subsatelllte dimensions (1.2 m radius In Table 1), the current Is 
determined by the lower electrode and It Is much smaller: 

Unreasonable balloon dimension would be needed to draw a current comparable 
to the Shuttle thermal currents. (This Is because the oxygen Ions are 
collected, at the lower electrode, at the Shuttle velocity v, whereas the 
electrons are collected, at the Shuttle, at their thermal velocity which Is 
nwch higher). 

An electron gun at the lower electrode, would allow reaching the current 
value of 0.19 amps determined by the Shuttle. To have more than that, again, 
an Ion gun at the Shuttle would be necessary. 

Similar considerations apply to configuration B (i.e., with the Shuttle 
below). In this case, referring to L * 100 km and always the Shuttle at 
h = 220 km, we find that 

^IO(s) 


for balloon radius of 


rjj -V 1 cm , 
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(as, this time. It Is the balloon that collects electrons). Thus, for 
reasonable balloon (or subsatelllte) dimensions (rjj -v* 1 m, for example). In 
a passive system and with an orbit Inclination reducing drastically the e.m. 
force, the current In the tether Is determined by the Shuttle collection and, 
hence. It Is given by 

1 -v 1.657 X 10"^ amps . (5) 

Notice that this same current, as written In Table 1, Is also obtained 
with all the available emf (2 * lO** volts) on the balloon for very small 
balloon dimensions. 

To Increase the current over the value (5) (taking advantage of the fact 
that for r^j > 1 cm, the balloon can drive higher currents then the Shuttle), 
one needs, for configuration B, an electron gun on the Shuttle. 

b) Limitation due to Internal tether resistance 


An upper limit for the tether current, which we call resistive current 
1|^, Is given by 


where aV Is the potential difference (due to V x b) between the electrodes and 
R the electrical resistance of the tether. Figure 1 gives curves of 1j^ versus 
wire radius r^, for steel and aluminum tethers. In the hypothesis of a 
naximum electromotive potential aV « VBL. Then, for exair^le, for a steel 
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tether of radius » 0.5 nm, we get 

Ip * 1.047 amps . (7) 

In comparing with the current values quoted above (valid for R « 0) we 
must however take Into account that, for Inclination of the orbit, such that 
aV becomes very small, the resistive current can also become very small. For 
example, (again for a steel tether of 0.5 rrm radius) the resistive current 
goes to the value 1.5 x lOr^amps, quoted above for a passive configuration A, 
for an angle e between and £ 

^ -v/ 0.82 . (8) 

Thus, for high Inclination orbit, one may be prevented, even using guns 
as explained before from reaching current values 1 'v 0.1 - 1.0 amps, useful 
for electrodynamic control. The only way out of this for cases where the 
V X B electromotive force becomes very small. Is to use guns which also 
accelerate particles, thus acting as voltage generators. Essentially, In this 
case, the aV of equation (6) Is supplied by the gun's acceleration with the 
result of avoiding the resistive limitations. 

c) External impedance 

The question of the external Impedance or, more generally, the coupling 
of the tether with the lower Ionosphere through magnetic flux tubes, has still 
to be considered quantatively (see Dobrowolny et al . . 1979, for preliminary 
conments). In this context one irajst recall the early work by Drell et al . . 1965, 
where the successive excitation of a series of transmission lines by the 
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motion of a large body In space was considered. We note only that the 
Impedance of the "Alfven wings" In Orel! et a1 . (1965) Is much smaller than 
any reasonable tether resistance. 
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1.14 Conclusions 


The study of electrodynamlc stabilization has shown that electrodynamic 
forces can be used to damp pendular oscillations of the tether. The damping 
can be applied to the total oscillation, or to a component of the oscillation. 
The most useful approach appears to be to use electrodynamlc forces to control 
the out'Of-plane oscillation. The perturbing forces out-of-plane are small, 
but there Is a resonant build-up of the amplitude, and tension control Is 
only weakly effective on the out-of-plane component. Relatively small currents 
can control the out-of-plane oscillation. There Is less need to control In- 
plane oscillations even though the perturbing force Is larger because of the 
lack of resonance. When necessary, tension control can be used effectively 
for the In plane oscillation. 


REFERENCES 


1. M.D. Grossl. "Orbiting tether electrodynamic Interactions" Final Report 
on Contract NASA-25077, Smithsonian Astrophys. Obs., April 1979. 

M. Dobrowolny, O.A. Arnold, G. Colombo, and M.D. Grossl. "Mechanisms of 
electrodynamic Interactions between a tethered satellite system and the 
Ionosphere." Smithsonian Astrophys. Obs., Reports In Radio and Geoastronomy 
No. 6, August 1979. 

3. S.O. Orell, H.M. Foley, and M.A. Ruderman, 1965. "Drag and propulsion of 

large satellites In the Ionosphere: An Alfven propulsion engine In space." 
Journ. Geophys. Res. 70, pp. 3131-3145. 


30 . 


2. Retrieval 


2.1 Introduction to Retrieval 

One of the potential sources of unstable tether behavior is retrieval 
of the subsatellite. This operation can Induce or amplify pendular oscilla- 
tions of the tether. A study has been done to see how electrodynamic forces 
could be used to help stabilize the retrieval process. Some new control laws 
using the reel motor have been devised which can be used together with electro- 
dynamic forces to dmnp oscillations of the tether. 

2.2 Negative Damping of Out-of-Plane Oscillations 

The coupling between the radial and out-of-plane variables introduces an 
«• 

acceleration given by 

• 

•• y* * 

A m .2 — i 
't'r r ^ 

If r is positive (deployment), the coupling always gives positive damping since 

the acceleration and velocity are of opposite sign. If r is negative (retrieval), 

the damping is always negative except when ^ is zero. The so-called “yo-yo" 

approach reels in the tether at the ends of the swing while ♦ is small and lets 

out line while ^ is large and ♦ is small. This is, of course, a time consuming method 

of retrieval. Electrodynamic forces can provide a way of damping this out-of- 

plane oscillation that is so awkward to deal with using just the reel motor. 

• • • 

For the in-plane variable, the coupling term is -2(0+y)rcos^. Unless the 
• . 
magnitude of 0 is larger than the orbital angular velocity y the acceleration 

is always in the same direction so that either negative or positive damping 

is possible using the reel motor. 
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2.3 Control Laws for Retrieval 


In order to study electrodynamic control during retrieval, the general 
equations of motion given In Appendix A have been Included In the small com- 
puter program used for the heuristic study of electrodynamic stabilization. 

The variables Included are the In-plane angle (e), the out-of-plane angle U) 
and the radial variable (r). The radial variable can be controlled either by 

. a tension control law, or a rate control law. If the tension Is controlled, 

• • # 

the program Integrates six quantities, e, e, r, and r. The second 
derivatives of each variable are given by the equations of motion. If a rate 
control law Is used, only five quantities are Integrated. The sixth quantity 
f Is specified by the rate control law. A rate control law does not guarantee 
positive tension. Runs done with the tension control law, 

TENSION - 7W + 4W i - 4W 2ic 
0 0 0 

(eq. III-l of NASA TMX-73314) show the expected behavior, namely a slow stable 
retrieval. The fastest possible retrieval Is one where the tether Is retrieved 
at a constant angle and the retrieval rate Is controlled such that the 
corlolls forces exactly balance the restoring forces In the In-plane direction. 

The equation of motion for e Is 

•• • ♦ • • • 

F9/m ■ rcos49+2(e+T)(rcos^-r^s1n^) + Sry^cos^coseslne. 

•• •• • • 

Setting gives 

• • • 

Zyr^ + 3ry^ cos 9^ sin 9. ■ 0 
o r r 


or 


3 

7 


r Y cos 9^ sin 9. 


(9) 


This Is the rate control law for maintaining a constant retrieval angle 9^ 

In the absence of any perturbing forces. A run done with the computer program 
using this rate control law with an Initial angle 9^, shows an exponentially 
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decreasing value of r with 0 ■ and ♦ » 0. 

A rate control law can also be used to provide damping In the In-plane 
direction. For small angles, the equation of motion for e Is 

•t 

m r 0 + 3 m r Y^e ■ F. 

In order to have critical damping we would like to add a term b0 to the 
equation where b is given by 

b ■ 2 /(mr) (3mrY^) 

. ( 10 ) 

■ 2/3 mrv 

The Coriolis term In the full equation of motion for 0 Is 

• • • • • • 

2 m (0 +y) rcoS(^ % 2m (0+y) r 

for small values of The value of r required to maintain a constant retrieval 
angle Is given by equation (9). To this value of r can be added an Increment 

ari defined by the equation 

• • • • 

b 0 ■ 2m (0 +y) ari 

In this way the corlolls force can be adjusted to provide the necessary damping. 
The damping part of the rate control law Is 

2 m (0 +y) 


Substituting the value of b from equation (10) gives 


Ari 


• • 

_ r Y 9 

• '■ I ’ I — 

2m (0 +y) 


( 11 ) 


A run has been made starting with 0 « 0. and applying a rate control law 

• • • 

r ■ r + Ar 
0 
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The tether angle Increases from e • 0 to e ■ smoothly, approaching 
exponentially. The retrieval remained stable to the end of the run. The above 
runs Indicate that a rate control law can control the In-plane behavior In the 
absence of any out-of-plane oscillations or any perturbing forces. 

The rate control law In equation (9) does not take Into account drag or 
other forces on the tether. In an analytical computation, the control law 
could take Into account the computed forces. In actual practice the forces 
may not be well known. The control law can Include a "restoring" force to 
counteract deviations of the tether In-plane angle from the desired value e^. 
The restoring force should be the same order of magnitude as the restoring 
force, 

3 ry2 cos ^ cos 0 sin 9, 

In the general equation of motion for the In-plane angle e. A factor k can 

be Included to provide additional stiffness. A term ltz In the rate control 

law can be confuted from the equation 

• • • • 

2 (0+y) ar 2 cos ♦ ■ -3 r cos ♦ cos a sin e k {e-0|,) 

which gives 

. 2 • k(0"6p) 

ara ■ - T cos e sin e — r : ■ 

^ (en) 

Adding all the rate control terms we have a rate control law given by 

• • • « 

r ■ r^ ♦ ari + ara 

In order to maintain critical dancing, the damping part of the rate control 
Is Increased by the factor The restoring force from ara has the effect 

of speeding up the retrieval In the case where the subsatelllte Is above the 
Shuttle, and slowing down the retrieval for downward deployment. 
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2.4 Cases Run 


In runs done for downward deployment, the restoring term halts the 
retrieval a few hundred meters from the Shuttle. For the upward deployment 
case, drag can be counteracted longer, but the dynamic eventually becomes un- 
stable. For either case, alternative techniques imjst be used to control the 
subsatelllte close to the Shuttle. Reasonable values of current In the wire 
(a fraction of an ampere) appear to be capable of controlling the out-of- 
plane oscillations In an Inclined orbit. Special considerations apply to 
equatorial orbits or orbits passing near the magnetic pole. 

Details of the cases described above are given In Monthly Reports 6 and 
7 dated April and May 1980. Figure 6 shows a sample plot from Monthly Report 
No. 6. The Shuttle altitude Is 220 km, the orbital Inclination Is 45 degrees 
and the tether length Is 10 km deployed upward. The Initial In-plane angle 
Is .4 radians (lagging the Shuttle) and the out-of-plane angle Is .035 
radians. A rate control law with damping Is used for the In-plane angle, and 
electrodynamic damping Is used for the out-of-plane angle. There Is no 
atmospheric drag. Part (a) of the plot 1$ the In-plane angle, part (b) Is the 
out-of-plane angle, part (c) Is the radial distance, and part (d) Is the damp- 
ing current. The first colimm In each plot Is the time In seconds and the 
second column Is the variable being plotted. The behavior Is reasonably 
stable at 12,000 seconds when the distance Is about 5 meters. 

2.5 Conclusions 

It appears that electrodynamlc damping can be useful for controlling the 
out-of-plane oscillations during retrieval. The In-plane angle can be con- 
trolled using a rate control law which Includes daim>1ng and restoring terms. 
Differential drag between the Shuttle and the subsatelllte will eventually de- 
stabilize the retrieval process when the tether length becomes too short. 

Other techniques must be used to control the subsatelllte during the last stages 
of retrieval. 
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FIGURE CAPTIONS 


F1gur« 6(a). 
Figure 6(b). 
Figure 6(c). 
Figure 6(d). 


In-plane angle (radians) vs. time (sec.). 
Out-of-plane angle (radians) vs. time (sec.). 
Radial distance (cm) vs. time (sec.). 

Damping currrent (amps.) vs. time (sec.). 
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Figure 6a 
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Figure 6c. 
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Figure 6d. 


3. Launcher 


3.1 Introduction to Tether Launcher 

The objective of this study Is to develop techniques for using the tether 
to launch satellites from the Shuttle. Satellites could be launched either by 
simply releasing them at the end of the tether or releasing them after acceler- 
ation by a rocket. The principle problem that must be studied In the dynamics 
of the tether Itself. The techniques nwst be such as to avoid breakage of the 
tether, or loss of tension or dynamic Instabilities that could lead to loss of 
control of the tether. Acceleration of the payload by a rocket Induces various 
tether motions that must be carefully studied. 

3.2 Subsatelllte Thruster Model 

The SAO Skyhook software package has been modified to Include thrusters 
on the subsatelllte. We have adopted a model with three Independent thrusters 
In the In-plane, out-of-plane, and radial directions. The directions of the 
thrusters are shown In Figure 7. The Shuttle Is at the origin 0. The x axis 
Is parallel to the vector from the center of the earth to the Shuttle, the z 
axis Is normal to the orbital plane, and the y axis Is perpendicular to the 
x -2 plane In the direction of the orbital motion. The vector r goes from the 



Figure 7. Coordinate system for thruster or subsatelllte. 
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Shuttle to the subsatellite. The In-plane angle 1s 8 and the out-of-plane 

A 

angle Is The In-plane thruster Is In the direction of 9, the out-of -plane 

A A 

thruster Is In the direction and the radial thruster Is In the r direction. 

The Skyhook software computes the directions of the In-plane, out-of- 

plane, and radial vectors for the Shuttle and components of the vector r along 
these directions. The sign convention Is different from that used In Figure 7, 
.so that the In-plane direction Is In the minus y direction (opposed to the 
orbital velocity) and the out-of-plane direction Is In the minus z direction. 
Ths sign convention Is the same In the radial direction. The equations for 

the unit vectors giving the directions of the three thrusters are; 

r«— x + '^y + — z 
y r r 

0«-t-x+^y 

p p 



where : 

r » /x^ + y^ + z^ 

p s /x2 + y2 


The Skyhook software has been modified to read In accelerations to be 
applied In each of the three directions and the time pfiiod over which the 
accelerations are applied. Provision has been made for specifying the 
accelerations during two time periods so that the system can be accelerated 
and then decelerated. 


3.3 Wire Oscillations Induced by Subsatellite Thrusters 

The following two sections present analytic expressions for estimating 
the transverse and longitudinal wire velocities that would be Induced when the 
subsatellite is accelerated by thrusters. For the transverse velocity a 
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simplified model Is used where the wire mass Is represented by a lump halfway 
between the Shuttle and subsatelllte. In fact, the wire motion consists of 
a transverse wave travelling down the wire. The present analysis was done for 
the purpose of order of magnitude estimation of the dynamics of the wire. The 
analysis of transverse oscillations assumes the duration of the acceleration 
and deceleration Is short compared to the natural period for transverse oscll* 
latlons. Transverse motion of the subsatelllte Introduces corloHs effects 
which can Increase the wire tension. In the case of short tethers especially, 
the Coriolis effects may be many times larger than the equilibrium value of 
the tension. 

If the wire tension Is nearly constant during the acceleration and de- 
celeration of the subsatelllte, the amplitude of the transverse wire oscilla- 
tion Is approximately equal to the distance travelled by the subsatelllte. 

This may be the case when the tether Is long, the velocity Is low, or the tether 
does not have enough time to stretch In response to the corlolls forces. In 
cases where the tension Increases by a large factor during acceleration and 
then drops again at the end of the deceleration period, the wire may be given 
a large velocity which cannot be restrained by the remaining tension in the 
wire. If the transverse velocity acquired Is larger than the velocity for a 
large amplitude oscillation. Instability will result. 

The formulas derived In the next two sections have been compared with the 
results of numerical Integrations done with the Skyhook computer program. The 
details of the comparison are given In Monthly Report No. 10 dated July 1980. 

3.4 Transverse Wire Velocity Induced by Subsatelllte Thrusters 

Approximate analytic expressions have been derived for estimating the 
magnitude of the transverse wire velocity that will be Induced during the 
acceleration and deceleration of the subsatelllte. In the diagram below, the 
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wire 1s represented by the mass m halfWay between the Shuttle at the bottom 
and the maneuverable subsatelllte at the top. 



Figure 8. 




The force F accelerating the wire mass Is 
F - Tx / (i/2) * 2 T x/i 


where T Is the wire tension. The velocity Imparted to the mass Is 



dt 


( 12 ) 


The wire tension T for an In-plane launch using equation A32 with r « 0 and 
r * i Is given by 


T 


Sm^iy^ + m^ie 


^ + 2mjt0Y 


(13) 


where y is the orbital angular velocity, m^ Is the mass of the subsatelllte, 

9 Is v/i, and v Is the subsatelllte velocity. 

Let us assume that the subsatelllte Is subjected to a linear acceleration 
or deceleration lasting from ti to tz and that the Initial velocity at ti Is 
Vi. We can then derive an analytic expression for Av^ that can be used for 
the acceleration, coast or deceleration period. The velocity of the sub- 
satellite as a function of time Is 
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V ■ Vi + a(t-ti) ■ (vi- ati) + at s vi + at 


04) 


where 


vi s vi - ati 


05) 


The angular velocity needed for equation (13) Is 

Vi + at 


• V 

•■r* 


The position as a function of time Is 

X ■ Xi + / vdt » Xj + ^ (Vi + at) dt 
"ti % 

where 1s the position at time ti. Performing the Integration gives 

X * X + Xi t + X 2 t^ 

0 


where 


s xi - Vi t - I ati 


Xi - Vi 

X 2 ■ y a 

Substituting equation (16) Into equation (13) gives 

T • T + Ti t + T 2 t2 
0 

• 2 

when T^ s 3mj2v^ ^ ^ 

Ti = 2 ~ a Vi + 2m^ a Y 

m_ 

T2 = f a^ 

Multiplying equations (18) and (20) gives 

Tx • + Cl t + Cj t^ + Cs t’ + C 4 S'* 


(16) 


(17) 


(18) 

(19a) 

(19b) 

(19c) 

( 20 ) 

( 21 a) 

( 21 b) 

(21c) 

( 22 ) 
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where 


T X 

(23a) 

0 0 

— • 
o 

X 

>-• 

—i 

X 

o 

(23b) 

T^ Xz + Ti Xi + T2 X^ 

(23c) 

Ti X, + T2 Xi 

(23d) 

T2 Xz 

(23e) 


Integrating equation (22) from ti to ta gives 

/* C C 

Tx dt • Cjt 2 -ti) + (ta-t?) + 

C, 4 4 C 5 5 

"5" (t2*tl) + "5” 


Substituting (24) and (12) gives the desired expression for the velocity 
Increment av^. 

In the special case where ti ■ vi ■ Xi ■ 0., the Integration of the 
product Tx can be done easily. The velocity, position, and angular velocity 
as a function of time are 


V ■ a t, 

X ■ ^ a t^ 


and 


0 ■ v/A ■ at/A . 

The tension given by equation (13) becomes 

T ■ Sm^Ay^ + m^A + 2m^A (^) y 

• 9 

■ 3m^Ay^ + m^a^t^ / A+2m^aty 


Multiplying by x gives. 
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Integrating with respect to time we have 

1 


Txdt ■ j m^atY^t^ + m^a^tV* + J- m^a^Yt** 


Substituting this Into equation (1) gives 


aVw ■ jj^ (aY^t^ + y + J- a^yt‘*/i 


(25) 


Since Av In equation (12) Is proportional to the Integral of Tx, It Is 
w 

desirable to minimize the amount of time the subsatellite spends at high velocity 
In order to reduce the excitation of transverse wire oscillations. 


3 . 5 Radial Velocity Induced by Transverse Thrusters 

Transverse velocity of the subsatellite produces a radial acceleration 
because of corlolls forces. From equation A32, we have for ♦ ■ 0 and r * z, 

r ■ Fr/m + ze^ + 2Z0 y + 3Zy^ (26) 

If the system Is close to equilibrium the first and last terms In (26) nearly 
cancel each other so that 

r ^ Z0^ + 2 z0y (27) 

If the velocity Is given by equation (14) we have, substituting equation (16) 
Into equation (27), 

r - + Rit + R 2 t^ (28) 

where 

(29a) 
(29b) 


R^ ■ ViVt + 2Viy 
Ri ■ 2 a(Vi/z + y) 
R2 ■ 
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(29c) 



Integrating equation (28) Me have 

tz 


AV. 


r dt 
tl 

^ (t|-ti) + ^ (t|-ti) 


( 30 ) 


The radial velocity av^ can be minimized by having the subsatelllte spend as 
little time as possible at high velocity. 

The analysis In this section can be applied to an out-of-plane launch by 
using the equation 

• 4 

T ■ Sm^iY^ + 

• • 

Instead of equation (13). Note that there Is no cross term In an out-of- 
plane launch. 


3.6 Tension Variations Caused bv Transverse Wire Oscillations 


The acceleration of the subsatelllte by thruster Induces transverse wire 
oscillations. These wire oscillations cause tension variations at the natural 
frequency for transverse oscillations by causing an acceleration r In the 


radial direction. The corresponding tension variation Is mr where m Is the 
mass of the subsatelllte. In Figure 9 below, the wire mass Is represented by 
a lumped mass halfway between the Shuttle and the subsatelllte. 



Figure 9. 
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The distance r between the ends of the wire Is 

r ■ 2 - y cos e) ■ z - £ cos 9 

The radial velocity Is 

• • 
r ■ i sin 9 9, 

And the radial acceleration Is 

• •• 

r ■ £(cos 9 92 + sin 9 9) . 

If we assume 


9*9 sin (I) t 
0 

from which we have 

9*9 (j) COS u t , 

0 

and 

9 ■ -9 sin u t • 

0 

then 

r ■ i(cose 92 cos^ mt - sin 9 9 ^ sin wt) 

0 0 

■ £ 9 (9 COS 9 COS^ u>t - sin 9 sIn ut) . 

0 0 

For small values of 9^, the factor In parentheses can be simplified. 

Making the substitution e^a^slnt^twe have 
•• 

r ■ £9 (I»^ (9 cos [9 sinwt] COS^ut 
0 0 0 

- sin [9^s1nwt] sin ut) 
y £9 (9 cos^wt - e^sln^wt) 

V 0 0 0 

Using the Identity cos2ut ■ cos^ut-sln^ut we have the approximate equation 
•• 

r £ 9^ cos 2 ut 

^ 0 


From the equation of motion for transverse oscillations. 
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nn^ y 0 ■ -2Te » 

we see that natural frequency Is u - where Is the wire mass and T 

Is the wire tension. 

3.7 Tension Variations Due to Longitudinal Oscmatlons 

The dominant source of tension variations In generally longitudinal 
oscillations of the subsatelllte at the end of the tether. The natural 
frequency u of the oscillations Is 

« - ^ 

where m Is the mass of the subsatelllte and k Is the stiffness of the tether. 

The stiffness k Is 

k ■ EA/i 

where E Is the elasticity of the tether material, A Is the cross sectional 
area and i Is the tether length. Neglecting tether mass, the oscillations occur 
about the center of mass of the Shuttle plus subsatelllte, so that one can use 
for i the distance from the subsatelllte to the center of mass. 

3.8 Transverse and Longitudinal Propagation Velocities 

Transverse acceleration of the subsatelllte at the end of the tether 
generates a transverse wave which travels down the wire at a speed determined 
by the tether tension and mass per unit length. The propagation speed Vt is 

where T Is the tension, p Is the density of the tether material, and A the 
cross sectional area. When there are substantial corlolls forces due to 
acceleration of the subsatelllte, the tension may be much greater than the 
equilibrium value. After deceleration or payload release, the tension may be 
very low. 
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The velocity for propegatlon of stress along the tether Is 
Vj^ ■ /ETp 

where E 1? the elasticity of the tether material and p Is the density. For 
kevlar with density 1.44 g/cc and Young's modulus 7 x 10" dynes/an^ the speed of 
sound Is about 7 x 10^ cm/sec. There will therefore be a delay In propagating 
tension changes from one end of the tether to the other. 

3.9 Coriolis Forces In-Plane and Out-of»Plane 

For small angles, the equilibrium Mire tension Is given by 
T ■ m (ra^ re^ ♦ Zr9y * 3nr^) 

There Is an asymmetry between 1n>plane and out*of*plane launches when the pay- 

• • 

load Is accelerated by thrusters. There Is a term 2rev which arises because 
the In-plane angular velocity Is parallel to the orbital angular velocity. 

There Is no corresponding term for the out-of-plane angular velocity. 

3.10 Coupling Between In-Plane and Out-of-Plane Oscillations 

Out-of-plane launches using thrusters to accelerate the payload produce a 

coi^nent of motion In the In-plane direction because of coupling between the 

variables. In-plane launches produce no out-of-plane motion as long as there 

is no Initial out-of-plane oscillation. The reasons for this behavior can be 

seen from the equations ?f motion. The coupling terms for the out-of-plane 
•• 

acceleration (^) are 

• • 

r sin a cos ♦ [e + y]^ 
and 

• • 

Era . 
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If 9 1$ 2 «ro, thtn the first term Is zero for all values of e because of 
the sin 4 factor. If 4 Is zero then 4 Is of course also zero so that the term 
2r4 causes no out-of>p1ane motion. The coupling terms for the In-plane 
acceleration (e) are 

2 (i+y) r cos 4 
and 

• • • 

-2 (e + y) r 4 sin 4 

The second term will produce an In-plane acceleration, whenever there Is 
an out-of-plane velocity except at the point where the out-of-plane angle Is 
zero. In addition If there Is a radial velocity, the first term will contribute 
to the In-plane acceleration. 

3.11 Reel Control Aloorlthtim for Pavload Release 

The al^rltNns described In this section are designed for releasing a 
payload at the end of a long tether In a stable manner. We assume that the pay- 
load has been deployed to the desired distance from the Shuttle and stabilized 
In a vertical configuration. The objective Is to release the payload In a 
manner that will not cause recoil of the remaining mass and loss of tension In 
the tether. A relatively simple approach Is to use the reel motor to Induce a 
longitudinal oscillation of the end mass, and release the payload when the ten- 
sion Is at Its minimum and close to the equilibrium value for the remaining mass. 

The initial study of this technique was done using a small program which 
Integrates the motion of two bodies connected by a massless tether. The action 
of the reel motor Is simulated by having the natural length of the elastic 
wire be a function of time according to the equation 
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Ola) 


Zq * U (t - t^), < t < 

^0 ■ ^c " ’ ^0^’ ^ ^ V 


(31b) 

(31c) 


The constant u Is chosen to be the same as the natural frequency for longitudinal 
oscillations of the end mass in order to maximize the induced oscillation. By 
choosing t^ - t^ equal to one half cycle, the effect is to pull the end mass 
toward the Shuttle during the first quarter cycle and then return the wire to 
its original length during the second quarter of a cycle before the end mass 
has a chance to return to its original position. 

Equation (31) has discontinuities in the wire velocity at t^ and t^ which 
tend to excite wire oscillations. An alternative to specifying length vs. time 
is to specify tension vs. time. We can use a tension control law of the form 


T ■ Tq + aT sin u (t - t^) t^ < t < t^ 

T * Tq + aT sin u (t^ ■ ^ ^ V 


(32a) 

(32b) 

(32c) 


When equation (31) is used, the minimum tension occurs at one half cycle. For 
equation 32, the minimum tension occurs at 3/4 of a cycle. Therefore, u in 
equation (32) has been set so that the period is 2/3 of that used in equation 
(31). There is considerably less excitation of wire oscillations using tension 
control. This should make ic possible to predict and control the tension at 
the time of release of a payload more accurately. 

More work needs to be done to develop algorithms for damping wire oscilla- 
tions during reeling operations and payload releases. 
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3.12 Other Methods for Avoiding Loss of Tension After Payload Release 


Section 3.11 describes a maneuver for avoiding loss of tension due to 
rebound after release of a heavy payload. In addition to the longitudinal 
oscillation of the payload after release there will be stress waves propagating 
along the wire. It should be possible to develop a control algorithm for the 
reel motor which will damp out these oscillations. When the end mass Is 
accelerated by thrusters, corloHs effects will tend to stretch the wire. 
Increasing the tension and causing energy to be stored In the wire which can 
cause rebound when the payload Is released and the remaining mass decelerated. 
The direction of the thrust could be controlled so as to provide a longitudinal 
component to cancel the corloHs forces. This longitudinal component could 
also be used to compress the wire to the point where the wire tension Is at 
the value anticipated after payload release. Deceleration after release will 
again reduce the equlllbrlian tension value. Th? thrust could be directed so 
as to bring the wire to an equilibrium condition without loss of tension. 
Additional work would be required to develop algorithms for controlling the 
thrusters so as to avoid loss of tension In the wire. 

Another possibility is to have a reel on the subsatellite that could 
deploy some wire to reduce the tension when a payload Is released. The amount 
of wire required would be on the order of the stretch In the wire which could 
be a kilometer or more. The shuttle could also release wire to reduce the 
tension. It would be necessary to anticipate the release because of the 
propagation delay along the wire. 

3.13 Long vs. Short Tether Launcher 

In a short tether, the equlllbrliwn tension due to the gradient of the 
gravitational and centrifugal forces Is small. If the end mass is accelerated 
In order to launch a payload the centrifugal forces generated may be large 
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compared to the equilibrium tension so that the tension goes through a 
large dynamic range. The angular velocity and angular displacement during 
the acceleration and deceleration phases may also be large. Therefore the 
velocities that can be achieved without generating Instabilities are rather 
limited. Only a small velocity is required to start the end mass circling 
around the Shuttle. 

In a long tether the equilibrium tension Is larger and the centrifugal 
forces are smaller so that the dynamic range of the tension Is smaller. The 
stretch of the wire Is larger so that more radial displacement can be accomo- 
dated without loss of tension In the tether. The angular displacement corre- 
sponding to a fixed lateral displacement decreases with tether length. The 
linear displacement and velocity of natural pendulum oscillations of the 
tether Is larger, so that larger subsatellite velocities are required to 
generate Instabilities. The mass of the tether itself Is larger compared to 
the mass of the subsatellite, and tends to be a stabilizing factor since more 
energy is required to excite motions of the tether. 

Considering all of the above factors, it would appear that long tethers 
are more effective and stable than short tethers. A short tether would 
probably have to circle the Shuttle to achieve large launch velocities. 

3.14 Cases Run 

Various cases have been run to study the use of the tether as a launcher. 
These have been done with and without thrusters, with short and long tethers, 
for In-plane and out-of-plane launches, with slack tethers, and with techniques 
for keeping the wire in tension. The details of the runs are described in 
Monthly Reports 8 through 13 dated June through November 1980. 

Some Initial runs were done with dancing between all of the mass points 
representing the wire. The behavior of such a system Is significantly differ- 
ent from a purely elastic wire because longitudinal oscillations of the wire 
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and end mass are suppressed. In practice, damping could be in^lemented at 
one or both ends of the wire. Therefore, most runs have been done with 
purely elastic forces or special control laws or damping at the ends. 

Because of the representation of the wire by discrete masses, the 
Skyhook program is not well suited to studying slack tethers. In a test run 
where the thrusters were fired in a direction such that the wire goes slack, 
tension spikes resulted when the individual mass points came back into tension. 
Increasing the nun^er of mass points might help give more reasonable results 
but it is difficult to conceive of any good way of modelling a slack tether. 
Effort has been concentrated on devising methodi of avoiding loss of tension 
in the wire. 

Runs have been done to verify the formulas and principles discussed in 
the previous sections. The expression for tension variations due to transverse 
wire oscillations is easily verified only when damping is used to suppress 
longitudinal oscillations. Without damping, longitudinal oscillations of the 
subsatellite at the end of the wire dominate the behavior of the tension. 
Launches using in-plane thrusters show no out-of-plane motion and larger 
Coriolis forces because of the addition of the tether and orbital angular 
velocities. Out-of-plane launches have lower coriolis forces and a slight in- 
plane motion developes because of coupling terms in the equations of motion. 

The excitation of longitudinal oscillations of the subsatellite and transverse 
wire oscillations is in reasonable agreement with the formulas developed. 

Runs using thrusters with short tethers show the large range of tension varia- 
l:ion and potential problems with unstable wire dynamics discussed in the section 
on long vs. short tethers. 

Test runs indicate that the release of a heavy payload, either in the 
static case, or after acceleration by thrusters, results in loss of wire ten- 
sion due to rebound unless special techniques are used to deal with the problem. 


Transverse acceleration of the payload Induces a radial oscillation because 
of corlolls forces. In one test run, release of the payload was delayed 
until the minimum of the tension cycle caused by the radial oscillation. 

This procedure avoided loss of tension, but larger transverse wire motion 
resulted because of the longer distance travelled by the end mass. 

In order to Investigate the behavior of the wire when thrusters are used 
some simulations have been done with five masses, three of which represent the 
wire. A two metric ton subsatellite at the end of a 200 km wire Is accelerated 
to 200 m/sec In 20 seconds and decelerated from 20 to 40 seconds. Figure 11a 
shows the behavior of the tension In each of the wire segments vs. time. Fig- 
ure 10 shows the numbering of the mass points. 



Figure 10, 


The tension plotted for each mass point Is the tension between that mass and 
the next higher numbered mass point. For the last mass, #5, It Is the tension 
to the Shuttle which Is mass #1. The Initial conditions are set up so that 
the tension Is equal between all the mass points. This is not an equilibrium 
condition since the tension Is greatest near the Shuttle and smallest near 
the subsatelllte. The long period variation of the tension is due to the 
longitudinal oscillation of the subsatelllte induced by the coriolis forces 
during acceleration and deceleration. The short period variation results 
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from longitudinal oscillations of the wire masses. The damping coefficient 
has been reduced to a negligible value In this run so that the longitudinal 
oscillations are not suppressed. There was no release of the subsatelllte 
and no loss of tension In the wire. 

Figure 11b shows the 1n-f'''ine motion of the wire. The In-plane axis Is 
expanded by about a factor of 50 relative to the radial axis. As the sub- 
satellite at the top Is accelerated It generates a transverse wave which 
travels down the wire to the Shuttle at the bottom. 

Runs have been done to study the release of a heavy paylaod without 
thrusters using reel control algorithms to avoid loss of tension In the tether. 
The simulations use 5 mass points, three of which represent the wire. Equili- 
brium initial conditions are set up for a 10 metric ton payload at the end of 
an 80 km wire. The length control algorithm described In Section 3.11 is used 
with t ■ 10 sec and t- » 195 seconds. 

0 t 

The tension as a function of time Is shown In Figure 12(a). The tension 
Is constant for the first ten seconds. If the speed of sound Is about 7 km/ 
sec In kevlar the propagation delay for stress along an 80 km tether Is about 
11.5 seconds. This seems to be In approximate agreement with the delay time 
seen In the plot. The maneuver sets up stress waves along the tether which 
cause tension fluctuations. After release of the payload, although the 
average tension Is near the desired value, wire oscillations cause loss of 
tension at various points along the wire. 

The potential tether materials being considered for the Skyhook system 
all have low hysteresis, so that there would be little damping of stress waves 
along the wire. However, since there Is strong compling between the action 
of the reel motor and the longitudinal oscillations. It should be possible to 
damp longitudinal oscillations by a suitable tension or rate control law at 
the end of the wire. A run has been done with the Skyhook program to test 
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the principle that longitudinal oscillations can be controlled by Introducing 
damping at the end of the wire. The run described In Figure 12(a) has been 
repeated with damping introduced between the subsatelllte and the wire mass 
adjacent to It. There was no damping In the other three wire segments, and 
equation (31) was used to control the natural length of the wire segment next 
to the Shuttle. The tension as a function of time Is shown In Figure 12(b). 

The wire oscillations Induced by equation (31) have been damped out, and the 
oscillations induced by the payload release at 195 seconds are also reduced 
with time. There is no loss of tension. 

In cases where the tether length and mass are larger and the payload is 
smaller, the effect of wire dynamics Is Increased. A run has been done using 
equation (31) for a 200 km tether with a two ton payload. The tension vs. 
time Is shown In Figure 13(a). The tension Is constant for the first ten 
seconds. The propagation delay should be about 28.7 seconds for stress waves 
along the tether. The tension variations along the tether are of larger 

relative amplitude and longer period than in Figure 12(a). 

Equation (31) has discontinuities In the v/ire velocity at t^ and t^ which 

tend to excite wire oscillations. An alternative to specifying length vs. time 
Is to specify tension vs. time. A run has been done using the tension control 
law given In Section 3.11. The tension vs. time Is shown In Figure 13(b). 

There Is considerably less excitation of wire oscillations than In Figure 13(a). 
This should make It possible to predict and control the tension at the time of 
release of a payload more accurately. 

More work needs to be done to develop algorithms for damping wire oscilla* 
tions during reeling operations and payload releases. 
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3.15 Conclusions 


The tether appears to be a useful device for transferring payloads be- 
tween orbits In an economical manner. In the case of release of a payload 
without acceleration by thrusters, the tether and the remaining mass at the 
end recoil resulting In loss of tension. This tension loss can be prevented 
and longitudinal wire oscillations can be damped using the tether reel motor. 
In the case of acceleration of the end mass by thrusters, both longitudinal 
and transverse wire motions are set up. Some techniques for avoiding loss 
of tension have been tested and other possible techniques have been discussed. 
Further work needs to be done to develop appropriate algorithms. Long tethers 
appear to be more stable than short tethers. As a result of the elasticity 
of the tether, the subsatelllte Is almost semi -free to maneuver at the end of 
a long tether without loss of tension. 
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FIGURE CAPTIONS 


Figure 11(a). 
Figure 11(b). 
Figure 12(a). 

Figure 12(b). 
Figure 13(a). 

Figure 13(b). 


Tension vs. time for each wire segment when the payload 1s 
accelerated by a thruster. 

In-plane wire configuration vs. time when the payload Is 
accelerated by a thruster. 

Tension vs. time In each wire segment using the length control 
algorithm. The tether Is 80 km long and the subsatelllte Is 
10 tons. 

Tension vs. time In each wire segment using the length control 
algorithm with damping In the wire segment next to the subsatelllte. 

Tension vs. time In each wire segment using the length control 
algorltNn. The tether Is 200 km long and the subsatelllte Is 
two tons. 

Tension vs. tine In each wire segment using the tension control 
algorithm. 
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Figure 11a 



Figure 11b. 
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Figure 12b. 
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Appendix A 

General Equations of Motion 

In this section, the general equations of motion are derived for a tether 
system modelled as a mass m at the end of a massless straight tether which applies 
a tension force on the mass m In the direction of the Shuttle. The Shuttle Is 
assumed to move at constant angular velocity In the x-y plane as shown in Figure 
Al. 



The position of the Shuttle is given by the vector p. In cases where the mass 
m is not negligible compared to the mass of the Shuttle, the Shuttle does not 
move with constant angular velocity. In this case, the vector p can be used 
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to define the position of the center of the system which does move at constant 
angular velocity. The position of the mass m relative to the center of the 
system is given by the vector r. The position P of mass m relative to the 
center of the earth is 

P = p + r (Al) 


The equation of motion of mass m is 

*• 

F = m P (A2) 

where F is the total force on the system. The force F is given by 

f - Fg + Fg (A3) 

where F^ is the gravitational force and F^ is all other forces such as the 
tension in the tether, atmospheric drag, etc. Differentiating equation (Al) 
twice gives 


P * D + r (A4) 

Since o and y are both constant, ^ is given by the usual expression for the 
centrifugal acceleration, and we can write 

p = -PY^ p (A5) 


where p is a unit vector in the direction of p. The second part of equation (A4) 
can be evaluated as 


rsiL(rr)= — (rr + rr) 

dt^ dt 

• ** 

=rr+2rr+rr 


(A6) 


In Figure Al, the unit vectors r, e, and ♦ form an orthogonal coordinate 
system. The derivative of any of these unit vectors is perpendicular to that 
unit vector and must therefore be in the plane defined by the other two unit 
vectors. This fact can be used to simplify the evaluation of equation (A6). 
From Figure Al we have the following expressions for the unit vectors. 
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A A A 

r * cos<>cosep + cos<>sin9Y + sin^k (A7a) 



9 » -s1n9o + C0S9y 

A A A A 

^ • -s1n(^cos9p - s1n(^s1n9y -*■ coS(^k 

(A7b) 

(A7c) 

# 

A 

We need to differentiate equation (A7a) to obtain an expression for r 
(A6). In doing the differentiation we use the relations 

In equation 


• 

- . ^ 
P = Y Y 

(A8a) 

and 

• 

^ * * 
Y = *YP 

(ASb) 


k » 0 

(A8c) 

to obtain the result 

r = -s1n^cos9<^p - coS(>sin99p + cos<>cos9yy 



-Sln<>Sin9<^p + COSKOS99Y - COS^Sin9YP 
+cos4i^k 

(A9) 

Equation (A9) must be a linear combination of equations (A7b) and (A7c). By 
inspection we see that the terms in each column of equation (A9) are multiples 

A A 

of either 9 or Therefore we have 


A ^ A -.A ^ A 

r = (>(^ + cos4»99 + cos4»y0 



= i(j> + COS4 i(9+y)0 

(AlO) 

Differentiating equation (AlO) to obtain r needed in equation (A6) we 

neve 


•• • 

A ••A ^A m m 9 ^ ••A 

r a (()({• + (jxf -sin<j)^(e+Y)e + cos<>90 

• 

+ COS<>(0+y)0 

(All) 


• • 

A A 

To complete the derivation of equation (All) we need e and Differentiating 
equation (A7b) with the help of equation (A8), we have 
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8 ■ - COS00O - Sln0YY 

- Sln00Y - COS0YP 
» - COS0(0-»4)p - Sln0(0+Y)Y 
■ - (0+y)(cos0p + sin0Y) (A12) 

Equation (A12) must be a linear combination of equations (A7a) and (A7c). Since 

A 

equation (A12) does not contain k, the only possible con^lnatlon Is to multiply 

A 

(A7a) by cos^, nwitiply (A7c) by sln^ and subtract to eliminate k. Doing this 
we have 

A A 

cos^r • s1n(^4 > 

cos^<>cos0p + cos^(^s1n0Y + cos«^s1n(^k 

_ ^ A A 

+ s1n2«>cos0p + s1n^^s1n0Y + s1n*cos4>k 

= COS0P + s1n0Y (A13) 

Substituting equation (A13) Into equation (A12) gives 

e 

A ^ ^ A A 

0 3. (6+y) (cos^r - sin^i^) 

A 

Next we evaluate which Is needed In equation (All). 

(A7c) with the help of equation (A8) gives 

^ ^ A A A 

♦ = - cosKos0(>p + s1n(^s1n00p 

- cos«frs1n0+Y - s1n^cos000Y + s1n(^sin0YP 

- s1n*ik (A15) 

Equation (A15) must be a linear combination of equations (A7a) and (A7c). By 
Inspection we see that each column of equation (A151 Is a multiple of either 

A A 

r or 0. Therefore equation (A15) can be written 

♦ * - K - s1n^00 - s1n^Y0 
= - *r - s1n(fr(0+Y)0 (A16) 


(A14) 

Differentiating equation 
- s1n(jicos0rr 
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Using equation (A14) and (A16) to complete the evaluation of equation (All) 
we have 

r * ♦ [-ir - sin^(9+Y)e] - sin(^i(^+Y)e 

+ cos^ee - cos<>(0+Y)2(cos^r - sin^<>) 

» r [-42 - 

+ e [-2isin^(e+Y) + ecos^] 

+ ♦ [♦ + coS(>sin4i{9+Y)^] (A17) 


Equation (AlO) and (A17) are the results needed to evaluate equation (A16) 

for r. Making the substitutions, we have 

•» 

•• A 

r a r r 

+ 2f [(|i<(i + coS(^(0+Y)e] 

+ r |r [-^2 - cos^<(i(9+y)^] 

+ 0 [-2|sin«>(0+Y) + 6cosi>] 

+ ♦ [♦ + coS(|»sin(|»(9+Y)^]| 

» r [r - r^^ - rcos2«^(9+y)^] 

+ 9 [2rcos<(»(9+Y) -2r|sin^(9+Y) + 0 rcos<j)] 

+ ^ [2r^ + r^ + rcos^sin^(0+Y)^] 

Equation (A5) gives p in terms of p. We must express p in terms of r, 0 , 

A 

and (^. The result can then be substituted along with equation (A18) into equation 
•• 

(A4) for Equations (A7a, b, and c) form a set of simultaneous equations that 

A 

can be solved for p. Part of the work has already been done by eliminating 

A A 

k from equations (A7a, and c) to obtain equation (A13) which contains only p and 

A 

Y. Jfeltiplying equation (A13) by cose and equation (A7b) by sine gives 


and 


AAA A 

cose(cos^r - sinc^^) = cos^ep + cosesiney 

A A A 

sinee = -sin^ep + sinecosey 


A-5. 



which can be subtracted to yield 

A A ^ * 

cosecos(>r - cos0s1n^(> - sinee * p (A19) 

Putting equation (A19) Into equation (A5) gives 

f* 

p * - rpY^ COS0COS* 

+ 0 PY^ s1n0 

+ cos0s1n^ (A20) 

Adding equations (A18) and (A20) gives for equation (A4) 

P » r[r-ri^-rcos2<>(0+Y)^-PY^cosecos*] 

+ 0[ercos^+2(0+Y)(rcoS(>-r|s1n^) + pY^sIne] 

+ ♦[r^+2ri+rcos*s1n4i(0+Y)^ + PY^cos0s1n<^] (A21) 

Substituting equation (A21) Into equation (A2) completes the derivation of 
the equations of motion for the variables r, 0, and ♦. All forces K must be 

AAA A 

resolved along the unit vectors r, 0 , and ♦. The tension force Is In the -r 
direction and needs no further discussion. Atmospheric drag and other forces 
such as radiation pressure are not In a fixed direction and must be evaluated at 
each Integration point. Equation (A7) gives the directions of the unit vectors 

AAA A A A 

r, 0, and as a function of p, y. and k. From Figure A1 we see that 

p = COSY? + sInY j (A22a) 

and 

Y • -slnY? + cosyJ (A22b) 

AAA 

Using equations (A7) and (A22), forces given In the 1, j, k coordinate system 

A A A 

can be resolved along r, 0 , and 

Neglecting higher harmonics, the gravitational force Is always In the -P 
direction. Therefore It Is useful to express the vector P In terms of r, 0 , 

A 

and Substituting equation (A19) Into equation (Al) gives 
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P ■ r{r + pcosecos^) 

A 

- 9ps1ne 
• ^pcosesln^ 

The gravitational force Is 

g p3 

where P Is the rwgnitude of 1^. From equation (A23) we have 

p2» r^+2rpcos8cos^ + p^cos^ecos^^ 

+ p^sln^e + p2cos20s1n2(j) 
a r2+2rpcos0cos^ + p^sln^e + p^cos^e 
or 

p3 . (r2+2rpcos0cos(^ + 


(A23) 


(A24) 


(A25) 


Substituting (A23) Into (A24) gives 

^ fiMn “ 

P _ wmi [r{r+pcos 0 cos(j>) 

3 p3 

- 0ps1n0 

- (^pcos8s1n(^] (A26) 

Combining equations (A26), (A20), and (A3) into equation (A2) gives 

F - .... 

mr[r-r42-rcos2*(0+Y)^-PY^cos0cos4 + GMP*^(r+pcos0coS(>)] 

+ m0[0rcos^+2(0+Y){fcos4-r^s1n4)+pY^s1ne-GMP’3ps1n0] (A27) 

+ mi>[ri>+2f^+rcos^s1n<^(l+Y)^+PY^cos0s1n(>-GMP“3pcos0sin4>] 
with P given by equation (A25). 

In cases where the tether length r Is small compared to the radius of the 
orbit p, the exact equations of motion given In (A27) can be simplified by 
expanding the term P’^. From equation (A25) we have 
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GMP"3 * GM(p2 + 2rpcos0cos(> ♦ 

■ 6Mp"3(l + 2 ~ cosecos^ + 

p p^ 

»6Mp-3(l - 3 j cosecoS(>) (A28) 

Using the relation GMp“^ ■ for a circular orbit, equation (A21) be- 
comes 

GMP"^ * - 3 - Y^cosecos^ (A29) 

P 

Substituting (A29) Into terms containing Gf1P‘^ In (A27) and neglecting terms 
containing r/p, we have the relations 

GM?“^(r + pcosecos^) rr^ + Y^pcosecos^ - 3rY^cos26cos2^ (A30a) 

Gf^P'^psIne pY^sIne - SrY^cosecos^sIne (A30b) 

6MP“^pcos0s1n^ Y^pcosesln^ - Sry^cos^ocos^sln^ (A30c) 

Substituting equation (A30) Into (A27) gives the approximate equations of motion 



mr[r-ri^ - rcos^(>(e+Y)2 + 1^2. 3rY^cos20cos^ifr] 

+ m0[0rcos^ + 2(e+Y)(rcoS(>-ris1n^) + 3rY^cos0cos<>s1n0] (A31) 

+ m^[r^ + 2ri + rcos*s1n*(e+Y)^ + Sry^cos^ecos^sln*] 

In equation (A31), p does not appear explicitly. It Is used In calculating 
Y, and It Is assumed that r, 0, and ^ are measured from a point at distance p 
from the center of the earth moving at constant angular velocity y. It Is useful 
to derive the equations of motion for small values of 0 and (^. The result Is 



mr[r - r^^ - ro^ -2r0Y - r-f^ + ry^ -3ry^] 
* m0[re + 2(e+Y)(r-r|^) + Sry^e] 

+ m^[r^ + 2ri +r^(J^+20Y+Y^) + 3r*Y^^] 


A-8. 



■ mr[r - r - 2r9Y - 3ry^] 

+ meCre + 2(S+Y)(r-m) + Srr^e] 

+ m^[r(> + 2f^ + r^(9^+29Y) + 4ry^^] (A32) 

If 9, I and f are also small, we obtain the simplest form of the equation namely, 

Fg Z mr(r - 3rf^) 

+ me (re + 3rr^9) 

+ ma(r^ + 4r*Y^^) (A33) 

From equation (A33) we see that the equilibrium tension Is 3mrY^, the 
natural frequency for In-plane oscillations Is y and the natural frequency for 
out-of-plane oscillations Is 2y. 
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Appendix B 

Orbital Coordinate System 



Figure B1 

In Figure A1, the x and y axes are taken to be In the orbital plane, and 

A A 

the z-ax1s Is the normal to the orbital plane. The unit vectors o and y are 
given by equations (A22) In terms of the unit vectors along the axes. In 
Figure Bl, the x'-y’ axes are in the equatorial plane of the earth with the x' 
axis taken along the ascending node of the orbit. The two coordinate systems 
are related by a rotation about fhe x* axis through 1 where 1 Is the Inclination 

A A A A 

of the orbit. The unit vectors 1, j and k are related to the unit vectors 1', 

A A 

j', and k* by the equations 


A 

r 

(81a) 

cos 1 J' ♦ sin 1 k' 

(Bib) 

-sin 1 J' ■►cos 1 ic* 

(Blc) 
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Substituting equations (Bla and b) Into equations (A22) gives 


A 

p ■ COS Y 

1* ♦ sin Y (cos 1j'+ 

Sin 1 k') 


■ COS Y 

r sin Y cos 1 j' + 

sin Y sin 1 k' 

(B2>) 

Y ■ - Sin 

A A 

Y r COS Y (cos 1 j' 

' + sin 1 k') 


■ - sin 

A A 

Y r ♦ COS Y COS 1 j ' 

cos Y sin 1 k‘ 

(B2b) 

Equations (B2a, 62b and Blc) define the unit vectors 

used In equation 

(A7) for 


an inclined orbit. 
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Append U C 
Initial Conditions 




Tht velocity of a mass point Is 


• wR^ (C2) 

The components of the velocity are 

i, ■ - V, xi / R, 

y, • V, X, / 

2^ ■ 0 (C3) 


In order for the system to be In equilibrium In circular orbit, the angular 
velocity u must be such that there Is no radial acceleration of the system. 

One condition of equilibrium Is that the sum of the gravitational and centrifugal 
forces on the total system must be zero. In addition, the sun of the tension, 
centrifugal and gravitational forces on each mass must be zero. The first condi- 


tions Is given by the equation 




(C4) 


1 

where r^ Is the distance of each mass m^ from the center of the earth, and w Is 
orbital angular velocity. Given the positions and masses of each point, the 
equilibrium angular velocity can be computed. The point In the system which Is 
at a radial distance p defined by the equation 


GH 


p 


2 


pw 


2 


which can be solved to give 



(C5) 


experiences no net radial acceleration. Points above and below o experience a 
net acceleration away from the center. 



The condition for equilibrium of the tension forces can be computed using 

•« »• •• 

the expression for r given in equation (A27). Setting r = 9 * ♦ * 0., and 
using the exact expression for the gravitational force since the tether length 
is appreciable, we have 


Fg/m = GM (pcos0cos^+r) P"^ - Y^pcosecos* - ry^cos2(^ (C6a) 


where 


P = + 2prcos9cos«> + r^ 


(C6b) 


and r is the distance from the center of the system. 

Using equation (C6), the force on each mass in the Skyhook system can be 
computed. The tension in each wire segment is equal and opposite to the sum of 
the forces on all masses which the segment supports. Starting at the end of the 
system, one can successively compute the tension in each wire segment. Once this 
tension is known, the natural length of the wire segment can be computed from 
the equation 

» Ar / (1 + (C7) 


where Ar is the distance between the ntasses, E is the elastic modulus, and A the 
cross sectional area jf the wire. Since the mass of each wire segment depends on 
the natural length which is a function of the tension, the solution for the 
equilibrium must be iterated to obtain a precise equilibrium configuration. 
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